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AN EXPERIMENTAL TEST OF VISUAL 
COMPARISON TECHNIQUE IN ESTIMAT- 
ING TWO DIMENSIONAL SPHERICITY 
AND ROUNDNESS OF QUARTZ GRAINS 
MELVIN A. ROSENFELD AND JOHN C. GRIFFITHS 


ABSTRACT. This experiment attempts to evaluate the visual comparison 
technique for the estimation of sphericity and roundness of quartz grains 
in sediments. The experiment was designed to enable the data to be 
statistically analyzed by means of analysis of variance—a powerful statisti- 
cal device which has a wide field of application in geology. 

The results of the experiment indicate that valid estimates of sphericity 
and roundness may be obtained by visual comparison with standards; but 
unless the variation due to differences between operators is suitably 
evaluated, visually estimated differences in sphericity and roundness can- 
not be attributed solely to differences between grains but are affected to 
an unknown degree by operator variation. The sympathetic relationship 
between sphericity and roundness, for example, appears on the basis of 
this experiment to arise from psychological bias on the part of operators 
and not from an inherent relationship between these properties of quartz 
grains. 


INTRODUCTION 


N attempting to estimate quantitatively the variation in 

the properties of sediments it is essential to evaluate the 
limitations of technique. Certain fundamentals must be fulfilled 
before a technique may be adopted as a routine for estimating 
variation in some particular property, and in testing the 
technique, the more restricted the test the more restricted the 
conclusions to be drawn from it. Ordinarily an experiment 
which purports to compile a set of measured data is somewhat 
circumscribed in its implications unless an extensive background 
of similarly measured reference data is available to which to 
relate the new experimental determinations. Before attempting to 


compile the extensive background it is, however, necessary to 
choose a suitable, efficient technique which will lead to meaning- 
ful measurements. Hence painstakingly thorough evaluation of 
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a technique is essential, and if this aspect is neglected it will 
inexorably lead to very considerable wasted effort. 

In the present case we are interested in attempting to evaluate 
limitations in the estimation of sphericity and roundness. The 
most rapid technique for both purposes and the only practi- 
cable technique for routine estimation of roundness of large 
numbers of grains is by means of visual comparison charts 
(Krumbein, 194la, and Rittenhouse, 1943), and if this ap- 
proach could be safely adopted, a large amount of data on 
different sediments could be compiled in a reasonably short 
space of time. 

Initially it is clear both from laboratory experiments (Krum- 
bein, 1941b) and actual measurements of quartz grains in sedi- 
ments (Pettijohn, 1949, p. 46-55) that the magnitude of 
variation measured on the conventional scale of 0-1 is likely 
to be very much less in the case of sphericity than roundness. 
On this basis, therefore, it is necessary to use greater precision 
in the case of estimation of sphericity than roundness and it 
may well be that two different techniques will be finally neces- 
sary. From the many experiments quoted in the literature (see 
Pettijohn, 1949, for extensive bibliography) it is also clear, 
and this is very fortunate, that it is much simpler to measure 
sphericity than roundness. 

Again it is reasonably clear from Pettijohn’s graphic por- 
trayal (1949, fig. 26, p. 54) that while the trends of relation- 
ship between sphericity and roundness are similar in different 
sediments differences between techniques cause as much dis- 
placement in the position of the trends as differences between 
sediments. 

Finally, we have learned by way of extensive experience that 
differences between operators performing the same tasks may 
well interfere with, and sometimes exceed, differences between 
sediments and/or techniques. 

We may conclude, therefore, that in order to test a technique, 
an experiment should be reasonably comprehensive and, as 
the description of any technique in the literature is presumably 
published with the intention that other investigators should adopt 
it, differences between investigators should always be evaluated 
in testing a technique. It may not be out of place here to 
emphasize that reproducibility as frequently reported can be 
very misleading ; for example, many chemical analysts will claim 
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quite justifiably that they can achieve duplicate or replicate 
analyses within enviable margins of precision, yet differences 
between analysts will far exceed this precision error and lead 
to very embarrassing results if not realized and evaluated before 
the compilation of the necessary extensive background of refer- 
ence data (e.g., Fairbairn, 1951). 

The present experiment is one of a series in which various 
techniques and various measurable properties of sediments are 
being tested with a view to ultimate adoption of suitable analyti- 
cal procedures for compiling the necessary extensive background 
of reference data on the magnitude of variation in the properties 
of sedimentary rocks. 


PROCEDURE AND DESIGN OF EXPERIMENT 


The details of the procedure used in the estimation of varia- 
tion in sphericity and roundness in quartz grains from sedi- 
mentary rocks by visual comparison with standard reference 
charts are adequately described in the original articles (Krum- 
bein, 194la, and Rittenhouse, 1943). In the present case 21 
grains of quartz were randomly selected from a beach sand and 
mounted in semi-permanent mounts; the 21 grains encompassed 


a variation of from 0.71 to 0.93 in sphericity and from 0.10 
to 0.80 in roundness. 

The published charts (see literature cited) were photographi- 
cally reproduced and brought to a standard size of 4” x 5” 
overall. Seven operators, all graduate students in the Division 
of Mineralogy of the Pennsylvania State College, took part in 
the test’ which was carried out in 1948 under the immediate 
supervision of one of us (M.A.R.). The experiment was initially 
designed to evaluate differences between grains and differences 
between operators but, by arranging the replicate runs of each 
operator on different days, variation from day to day was 
also included. In terms of experimental design the differences 
between grains, differences between operators and differences 
between days are termed main effects, and subsidiary factors 
comprise three first order interactions (Youden, 1951), namely, 
operators and days, operators and grains, and grains and days. 
Finally, a second order interaction of operators and grains and 
days completes the list of effects included in the experiment. 


1The seven operators were C. T. Bressler, E. Cressman, D. B. Doan, 
R. L. Folk, M. A. Rosenfeld, A. Swineford, and C. E. Weaver. 
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marks are appended by way of brief explanation of the mean- 
ing of each source of variation; for greater details the interested 
reader may be referred to relevant statistical texts (Fisher, 
1948; Snedecor, 1946; Youden, 1951, etc.) and, for examples 
of similar procedures applied to geological problems, to a grow- 
ing number of publications (Swineford and Swineford, 1946; 
Chayes and Fairbairn, 1951; Miller, 1949; Griffiths and Rosen- 
feld, 1950; Rosenfeld and Griffiths, 1950). 

The main effects arising from differences between grains, 
operators and days respectively may appear to carry fairly 
obvious implications, but the use of analysis of variance is best 
understood when it is realized that each main effect is used to 
test a separate null hypothesis (Fisher, 1949). Thus we may 
set up a null hypothesis concerning grains as follows: 

There are no differences between grains in sphericity or round- 
ness; as a corollary the entire set of measurements is then a 
random sample of an infinite population of such measurements 


in which there are no differences in sphericity or roundness as the 
case may be. 


Our test consists of an attempt to disprove this hypothesis. 
The setting up of a null hypothesis actually implies that we 


have defined a population of which our set of measurements is 
a random sample, and our test comprises computation of statis- 
tics, the mean and variance, from the sample and estimation 
of how likely such statistical values are to arise on the basis 
of our hypothesis about the population. 

In the present case, as implied by the title “analysis of vari- 
ance,” we shall be testing the variance around the means and 
using the mean and variance but especially the latter as a basis 
of judgment of whether the sample of measurements has indeed 
been drawn from a population in which the differences in 
sphericity are non-existent. Each main effect leads to a separate 
null hypothesis and, of course, to a separate population for 
testing. The interactions are, however, somewhat more intricate ; 
we may compare, by way of illustration, the main effect dif- 
ferences between operators and the first order interaction 
operators times grains. There may be and, indeed, usually are 
differences between operators in which one operator always 
estimates, say, a high value of sphericity whereas another 
measures a low one (see, for example, operators in figure 5). 
In addition some operators estimate consistently, i.e., always 
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yield high or always low values, whereas others are inconsistent, 
sometimes measuring high and at other times measuring low 
values. The differences between operators is estimated as the 
main effect whereas the inconsistency of operators from grain 
to grain is measured by the interaction term operators times 
grains. It is reasonably simple to extend this argument to 
variations of operators from day to day and so to account 
for operators times grains and operators times days first order 
interactions. There appears to be no obvious physical meaning 
to the first order interaction of grains times days principally 
because it is difficult to conceive of a grain varying from day 
to day or being inconsistent in this sense; nevertheless, this 
term is mathematically separately evaluated here although it is 
customary to include the variance from this source with “error” 
in the usual case. The reason for this departure from conven- 
tion is twofold; firstly, it is simpler to explain the analysis in 
this way and, secondly, we are in such a dubious position 
geologically in terms of evaluating what is and what is not 
important that it behooves us to be extremely cautious in at- 
tempts to simplify by jumping to conclusions where we have 
no prior experience to guide us. Although, therefore, we can- 
not at present imagine a meaning for such a source of variation, 
it is as well to be aware of its possible eventuality in future 
experiments ! 

Finally, the second order interaction entitled grains times 
operators times days includes all variation from unassigned 
sources. This then becomes our basic error term for the experi- 
ment. It may be as well to add that this term is capable of 
physical interpretation in experiments possessing a more com- 
plex design and has been of value in evaluating certain agri- 
cultural experiments. In the present case, however, we do not 
attempt to interpret the physical meaning of the term but ac- 
cept it as a measure of experimental error. The magnitude of 
this error determines the sensitivity of the entire experiment. 


TREATMENT OF THE DATA 


The use of analysis of variance technique carries with it cer- 
tain implications concerning the variation in the measured 
variable (Eisenhart, 1947); ideally the frequency distribution 
of the measurements should follow the normal distribution but 
considerable departure from this ideal is permissible in practice. 
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Nevertheless, it is necessary to ensure that the data are suitable 
for this kind of statistical treatment, and certain standardized 
tests exist for this purpose. Furthermore, it is always advisable 
to be aware of the frequency distribution of the measured vari- 
able because this information is of considerable importance both 
in understanding the significance of the measurements and in 
evaluating the usefulness of the technique. 

The scales of measurement of both sphericity and roundness 
extend from 0 to 1 and are, therefore, “percentage” or propor- 
tion type scales, whereas the limits of the ideal normal distribu- 
tion extend from minus to plus infinity; as a result, the fre- 
quency distributions of sphericity and roundness may show 
cutoffs at zero and unity. The recommended practice is to 
transform the variable from percentage or proportion to angu- 
lar values (the inverse or arc sin transformation), and tables 
to facilitate this transformation may be found in many stand- 
ard texts (Fisher and Yates, 1948; Snedecor, 1946; etc.). 

In using analysis of variance techniques a common error 
variance is computed as the basis for comparison of all main 
and subsidiary effects, and for this reason it is necessary to 
ensure that the variance and especially the error variance is 
not related to the mean; the independence of mean and variance 
in the normal distribution is, of course, a fundamental prop- 
erty of such a frequency distribution. In order to ensure as 
near fulfillment of these prerequisites as possible, it is necessary 
to examine the frequency distribution of the variable in some 
detail, and should any doubt concerning the requirements exist, 
it is advisable to transform the variable to a scale in which 
as many as possible of the assumptions involved in using 
analysis of variance technique are fulfilled.” 

Sphericity measurements.—The 441 values of sphericity com- 
piled for this investigation are listed in table 2 and are graphi- 
cally displayed as a histogram in figure 1, together with the 
computed normal curve possessing the same mean and stand- 
ard deviation. The importance of the departures may be 
estimated in several ways and here the chi-square test was 

2A more detailed and formal discussion of this and other aspects of 
analysis of variance is included in a series of articles in Biometrics Bulletin 
(Eisenhart, 1947; Cochran, 1947; and Bartlett, 1947), and the use of the 
inverse sin transformation is more fully ventilated in “Selected techniques 


of statistical analysis,” Fisenhart, C., et al., McGraw-Hill Book Co., 1947, 
part 3, p. 395-416. 
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Fig. 1. Histogram of raw data sphericity values. 


used (Dixon and Massey, 1951, p. 190). The computed chi- 
square value yields a probability of less than .001; this may be 
interpreted to signify that if our 441 values were a random 
sample from a normal population of the same mean and 
standard deviation, such a chi-square would occur less than 
once in every 1000 samples. We judge, therefore, that our 441 
values depart significantly from a normal distribution. 

If we now calculate the mean and variance for each of the 
21 grains based on the 21 independent estimates for each grain, 
we find that the variance shows no obvious relationship to the 
mean. For our present purpose and based on this limited set 
of data we may deduce that, although the departure from 
normalcy is significant, the effect on the levels of significance 
in the analysis of variance and particularly the F test (Snede- 
cor, 1946) is not serious. Nevertheless, for purposes of explora- 
tion, each value was transformed to its are sin equivalent and 
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P< 


No. of Grains 


Arc sin Sphericity 


Fig. 2. Histogram of sphericity values—data transformed to arc sin. 


the analysis of variance performed on both raw and trans- 
formed data. The comparison of the results follows later 
(table 5). 

As a matter of interest the histogram of the transformed 
data is given in figure 2, together with the normal curve 
computed for the same mean and variance. The chi-square 
test yields a probability value of less than .001. It seems clear 
that the are sin transformation in this case did not suffice 
to normalize the data. It should be noted from the frequency 
distribution of the raw data that the mean value of 0.858, the 
almost symmetrical histogram, and the lack of cutoffs at either 
end suggest that no transformation is necessary in this case. 

Roundness measurements.—The roundness data based on the 
441 estimates is given in table 3 and plotted as a histogram 
in figure 3. This distribution tested against the normal curve 
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Fig. 3. Histogram of raw data roundness values. 


computed for the same mean and variance yields a chi-square 
of 41.84 and P <0.001, again indicating a non-normal distri- 
bution. Furthermore, a plot of the mean roundness against 
variance for each of the 21 grains yields a definite trend (fig. 
4). In this case the association is estimated by the correlation 
coefficient r = —0.4272 which, on the basis of 21 sets of 
values, yields a P of <0.1> 0.05 (Fisher, 1948, table 6, p. 
209). Although not significant at the conventional level, such 
a degree of association is clear indication that the mean and 
variance are not independent in this set of data. Transformation 
to inverse sin equivalents failed either to improve the normalcy 
or to render the mean independent of the variance. It was 
decided, therefore, to use Bartlett’s approach (1947, p. 39) 
and to calculate a correction from the data themselves. The 
trend of mean versus variance was assumed linear and com- 
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Fig. 4. Plot of mean roundness versus variance for 21 observations 
per grain. 


puted by least squares (see fig. 4) yielding s* = 0.01107- 
0.01208x, where s* = variance, and x = mean of the round- 
ness values for each grain respectively. Following Bartlett's 
recommendations (op. cit. p. 39) the linear relationship was 
integrated to yield the required transformation: 


Ps = V0.01107 — 0.0121P (1) 


where P = measured roundness value and Pa = adjusted round- 
ness value. Again the analysis of variance was computed on 
the basis of raw data, inverse sin and adjusted values as a 
basis for comparison (table 7). It should be noted that, using 
the computed formula given in equation (1) above, any round- 
ness value greater than 0.7 will yield a negative adjusted value, 
i.e., this transformation is not suitable for adoption other than 
with the data herein recorded. Needless to say, the regression 
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adjusted roundness values show a mean independent of the 
variance, 

In conclusion, it should perhaps be emphasized that the 
requirement of normalcy, while desirable, is not so essential as 
that the variance should be independent of the mean, and this 
is particularly true if, as in the present analysis, one is more 
interested in the estimation of components of variance than 


in the difference between means (see, in particular, Eisenhart, 
1947, p. 19). 


SOURCES OF VARIATION IN SPHERICITY 


The analysis of variance.—The analysis of variance of the 
sphericity data in terms of the transformed (arc sin) values is 
summarized in table 4. The steps in the analysis may be out- 
lined as follows: 

From the components of variance (see mathematical model 
table 1) it can be seen that each source of variance may be 
tested against its appropriate error term; thus, the first order 
interactions are tested against the second order interaction. 
The actual test is performed by comparing the mean squares 
and computing the variance ratio or F value (Snedecor, 1946, 
p. 219). For example, if we divide the mean square for the 
days times grains interaction by the mean square for the second 
order interaction we are evaluating the variance ratio, F=(s*,+ 
7s* nq)/s’e; returning to our initial argument, if there are 
no differences, then these two independent estimates are esti- 
mates of the same variance. In the ideal case, s* pg is zero and 
the variance ratio becomes (s*,/s*,) which is, of course, unity. 
Departures from this ideal value which arise from chance 
causes or from random sampling of a homogeneous population 
are estimated by the appropriate F distribution. In the present 
case the F value of 1.04 is a near estimate of unity; in fact, 
the 20% point of the F distribution with 40 and 200 degrees 
of freedom is approximately 1.208, and hence we would expect 
such an F value more than 20 times in every 100 of such 
tests in sampling from a homogeneous population. It would 
appear, therefore, that on the basis of the present evidence, 
variation from inconsistency of grains from day to day adds 
nothing more to the variation than that contained in the error 


8’ Computed using Fisher’s formula, Fisher and Yates, 1948, p. 34, foot- 
note to table. 
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Fig. 5. Variation in sphericity for each operator, and for each operator 
on each day. 


variance, In such a case we are justified in combining these 
two estimates to obtain a more efficient estimate of the error 
variance. The pooled estimate is found by adding the sums of 
squares for each term and dividing by the summed degrees of 
freedom; we obtain a mean square of 3.329 with 280 degrees 
of freedom and this value is recorded at the foot of table 4. 

The next step is to test the interaction operators times days 
and this yields an F value of 4.056 when tested against the 
pooled error. The 0.1% point of the F distribution for 12 and 
280 degrees of freedom respectively is 3.02, and hence we would 
expect such a large F value as that found very much less 
often than 1 time in 1000. We cannot but conclude that our 
initial hypothesis of no differences in the population is probably 
false in this case. Hence this interaction, the inconsistency of 
operators from day to day, contributes a significant amount 
of variance to the total. The variation of each operator for 
each day is illustrated in figure 5. Similarly the operators 
times grains inconsistency is significant at the 0.1% level. 

In attempting to evaluate the main effects we have a choice 
of two tests and each test supplies somewhat different informa- 
tion. Thus in the case of differences in sphericity from day to 
day the logical error term might be days times grains or days 
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times operators. Let us assume for the moment that differences 
between operators are omitted from consideration (e.g., the 
tests were performed by a single operator) ; then we may test 
differences between day means against days times grains or 
the more efficient estimate of the pooled error. The resultant 
F value of 3.88 exceeds the 5% F value of the equivalent F 
distribution (F,, = 2.99, F,,; = 4.60) but is less than the 
1% value. It seems clear that the differences between days 
significantly exceeds the experimental error and this suggests 
that it may be as well to include in any subsequent experiment 
a more extended period to evaluate differences between days 
before concluding that this source of variation may be neglected. 

The second test of differences between the means for days may 


be made by comparison against the inconsistency of operators 
from day to day, and this yields an F value of less than unity. 
In effect, then, the contribution to variation of the inconsistency 
of operators from day to day reduces the contribution due 
to differences between days to negligible proportions. 

The conclusion to be drawn from these two tests is that with 
a number of operators the inconsistency of operators from day 


to day is sufficiently large to overshadow differences due to 
measurements performed on different days; whereas when 
operator difference is removed, such as occurs when the tests 
are performed by a single operator, the day to day differences 
are sufficiently large to necessitate separate evaluation. These 
findings need confirmation from a more extended test over a 
longer period before final acceptance. 

Similarly, differences between operators may be tested both 
against operators times days and operators times grains inter- 
action. In both cases the results are substantially similar, both 
I values being significant beyond the 1% value of the cor- 
responding F distributions.* Differences between operators are 
clearly established. The variation in operator means is con- 
tained in the graph of figure 5 (solid line). 

Differences between grains tested against operators times 
grains interaction are also highly significant, the resultant F 
value far exceeding the 0.1% value of the F distribution for 
20 and 120 degrees of freedom. The differences between grain 
means are illustrated in figure 6. 


4In fact, operator differences tested against operator x grains inter- 
action is significant beyond the 0.1% level. 
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Fig. 6. Variation in mean sphericity and roundness among grains. 


Comparison of the analysis of raw and transformed data.— 
The analysis of variance performed on the raw data and on 
the data transformed to its arc sin equivalent is displayed in 
table 5 in terms of the resulting F values and their correspond- 
ing significance levels. It can be seen that only minor changes 
in levels of significance occur and no differences would arise 
in the conclusions drawn from the analysis of the data based 
on either set of values. This is only to be expected in this 
case where the original frequency distribution does not suffer 
from the kind of aberration which the inverse sin transformation 
is used to ameliorate. Nevertheless, it is comforting to know 
that no great change occurs in using the transformation where 
it is not essential; if, for example, two sets of data were to 
be compared, one of which necessitated the are sin transforma- 
tion and the other did not, such a property of the transforma- 
tion is a distinct advantage. 

Estimation of components of variance in sphericity.—One 
of the advantages of the use of analysis of variance in the 
treatment of data arises from the fact that a suitable estimate 
of the variance contributed by each source may be obtained 
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from the data. The items contributing to the variance of each 
factor together with their weighting coefficients are listed op- 
posite the respective sources in the mathematical model of 
table 1. It will be noted, for example, that the variance arising 
from variation of operators from day to day—operators times 
days interaction—is composed of a contribution from error 
(s*.) and a second contribution from the interacting items 
(s*yp). The second component in this variance is based on the 
21 grain means and the coefficient g is used to make the esti- 
mate conform with the error variance.® Thus we have an esti- 
mate of the operators times days interaction as s*, + gs*xp; 
we also have an estimate s*, for the error variance, and so the 
magnitude of the variance for the interaction (s*xp)) may be 
obtained by simple arithmetic. The variance for each term in the 
analysis is listed in table 6. 


It can be seen that the largest contribution to variance, or 
the greatest source of variation, arises from differences between 
grains, and this is almost twice the next largest contribution 
from unassigned sources—the error term. The third largest 
contribution is from differences between operators and this is 
about one-fifth of the variance from differences between grains. 


These components of variance clearly indicate that there 
must be many sources contributing to error, and from the 
magnitude of this residual variation there are other effects which 
could be profitably introduced into the analysis. In addition, 
this analysis of the variance components suggests that should 
the range of variation in sphericity between grains be reduced 
then the differences between operators will show a relative 
increase. 


Again it can be seen that the variance arising from incon- 
sistency of operators from grain to grain is about twice as 
large as that from the inconsistency of operator from day to 
day. The variance contributed by differences between days is 
negligible. 


‘This rule arises from the relationship between the variance of means 


and the variance of the raw items, ie., ns- — s where s— is the 
i x x 


variance between the means, s_ is the variance of the original items, 
x 


and n is the number of items used in calculating each mean (Crump, 1946). 
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The standard deviation derived from the variance is also 
listed in original units in the column headed “decoded s” but 
it must be remembered that the relationships between the 
standard deviations are not as simple as those between the 
variances. The coefficient of variation is computed from the 
error standard deviation (square root of the error mean square) 
and the grand mean of the entire experiment® and is considered 
to be a reflection of the experimental control; in general, a 
coefficient of variation less than 20% is reasonable and less 
than 10%, good. 


The comparison of variance components based on estimation 
using the raw data and the transformed data are also included 
in table 6. The relative magnitudes are certainly not the same, 
but the conclusions based upon them are, so that, in this case, 
the use of the transformation is evidently not necessary. 


It must be emphasized that this simple technique for the 
evaluation of the variance contributed by each factor in the 
analysis is based upon assumptions concerning the contributions 
to variance in the original data. One such assumption, implicit 
in the manipulations leading to separate evaluation of the 
variances associated with each factor, is the additivity of the 
variance from each source. If the original frequency distribu- 
tion were normal such conditions would hold, but in non-normal 
distributions it is necessary to demonstrate the independence 
of mean and variance for such conditions to hold. Transforma- 
tion of the variable is one approach to this problem and ac- 
counts for the attempt to use the inverse sin conversion in this 
investigation (Bartlett, 1947; Eisenhart, et al., 1947). In the 
case of sphericity determination, however, although the fre- 
quency distribution is non-normal, the mean does not appear 
to be related to the variance in any systematic manner. The use 
of the are sin transformation is not therefore mandatory in 
this case. In general, however, it is as well to ascertain empiri- 
cally that additivity is likely before the analysis of components 
of variance is attempted (Crump, 1946). 


E.M.S. 

® According to the formula Cv% = x x 100; where Cv= coefficient 

of variation, E.M.S.—error mean square, and Xg is the grand mean. 
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SOURCES OF VARIATION IN THE VISUAL ESTIMATION OF ROUNDNESS 


The analysis of variance.—The analysis of variance of round- 
ness in terms of the regression adjusted data is illustrated in 
table 7 and it will be unnecessary to repeat the steps in detail 
because the same procedure is followed as that described in the 
analysis of the sphericity data. The results of the tests and 
minor modifications appropriate to the different data are out- 
lined in the following paragraphs. 

Testing the days times grains and operators times grains 
mean square against that for error yield F values which are 
not significant. The variance for these terms may then be 
pooled, and the pooled error is recorded at the foot of table 9. 
Little is gained by this manipulation, and for subsequent tests 
the unmodified error term is used for testing the grains times 
operators interaction. This is significant beyond the 0.1% point 
of the F distribution for 120 and 240 degrees of freedom, and 
the F value would only be increased by using the pooled error 
term. Clearly the inconsistency of operators from grain to grain 
makes a considerable contribution and becomes the basic error 
for comparison with main effects. 

The difference between days yields an F value less than 
unity when tested against either the pooled error or the grain 
times operator interaction. Differences from day to day are 
therefore not significant. 

Differences between operators and between grains are each 
significant beyond the 0.1% level of the respective F distribu- 
tions. Hence these differences are highly important; the magni- 
tudes of variation for these two main effects are included in 
figures 7 and 6 respectively. 

Comparison of the analysis of raw and transformed data.— 
The comparison of the results of analysis of variance using 
raw data, are sin and regression adjusted values is summarized 
in table 8. No difference in the conclusions to be drawn from 
the analysis would arise whichever kind of data were followed. 


This serves to illustrate the degree of relaxation of the assump- 
tions of normality which may be permitted without affecting the 
levels of significance in using the analysis of variance procedure 
for testing the differences between means. 


Estimation of the components of variance in roundness.— 
An estimate of the contribution to variation arising from each 
term in the analysis is recorded in table 9, and the results 


(raw data) 


Roundness 


Mean 
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roundness 


Operator 


Fig. 7. Variation in mean roundness among operators. 


of this computation using the raw data, the are sin transformed 
data and the regression adjusted data are included for 
comparison. 

The qualitative treatment of these values by comparing their 
relative magnitudes in terms of each of the respective scales 
would hardly differ in any important respect, but it should be 
emphasized that for further manipulations the most reliable 
estimate is that based on the regression adjusted data. This 
conclusion is predicated upon the fact that in both raw and 
are sin data the mean is related to the variance and conditions 
upon which the computations of the estimates of variance are 
based do not hold. In the case of the regression adjusted data 
the mean is independent of the variance and additive relations 
between the components of variance are most likely. The sub- 
sequent discussion utilizes all three estimates, but attempts to 
compute the number of items necessary for the reduction of 
any variance component for future experiments must rely solely 
on the regression adjusted data. 
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The main source of variation arises from the differences in 
roundness between the different grains, and this term is at least 
twice that of the next largest source of variation, the error 
term. The variation associated with the error term is unassigned 
so that we may consider that the experiment has reduced the 
error variation to less than half of the variation we wish to 
determine. Evidently differences in roundness between grains 
may be established by the visual comparison technique. 

The third source of variation in order of magnitude arises 
from differences between operators, and this is approximately 
one-sixth of the variation arising from differences in roundness 
between grains. The next source, the interaction term, operators 
times grains, also contributes an important part of the total 
variation. In this experiment the basic error term for the com- 
parison of main effects is this interaction term, and hence the 
inconsistency of operator from grain to grain determines the 
level of sensitivity of the visual comparison technique. It can 
be realized from the form of the analysis that differences between 
operators can be segregated but their inconsistency determines 
the magnitude of the difference in roundness between grains 
which can be detected by this technique. It should be emphasized 
that any reduction in the variation in roundness between grains 
will result in a relative increase in magnitude of the incon- 
sistency; hence in any cooperative experiment this term should 
always be included and evaluated as a measure of sensitivity of 
the experiment. 

The other sources, differences between days, and the respec- 
tive interactions, operators times days, and days times grains 
do not contribute in any appreciable degree to the variation 
and may be included in the error term without inflating that 
item if future experiments confirm the present results. 

The recording of the decoded standard deviation in table 9 
is included merely as an appropriate guide to the relative im- 
portance of the several factors. In these experiments as shown 
by the coefficients of variation the use of the regression adjusted 
data has resulted in an improvement in the experimental control. 


DISCUSSION AND CONCLUSIONS 


It should be clear from the preceding analysis that valid 
differences in both sphericity and roundness of quartz grains 
may be established by the use of visual comparison techniques. 
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Fig. 8. Plot of mean sphericity versus mean roundness for 21 observations 
per grain. 


It is, furthermore, confirmed that variation in roundness is 
larger than that in sphericity. 

The magnitudes of the variation may be compared by means 
of figure 6, and the statistic “least significant difference” (LSD) 
may be considered as a very approximate measure of the least 
differences which can be detected in each case.’ It would appear 
to be prudent to measure sphericity rather than estimate it, 
both because measuring techniques available are simple and 
rapid and because the entire range of variation in sphericity 
is not likely to be much larger than that included in this 
experiment. 

In any future experiments some measure of operator varia- 
tion should be included, and in cooperative experiments includ- 
ing more than one operator the inconsistency of operators from 
grain to grain should also be evaluated. It must be emphasized 
that where this has not been carried out the differences in 
either sphericity or roundness may not be solely ascribable to 
geological factors but include experimental or technique error. 

7 This measure should not be taken too literally because it achieves precise 
meaning only under very restricted conditions and particularly where 


the frequency distributions of the measured variable are homogeneous and 
closely approximate normality. 
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Fig. 9. Plot of mean sphericity versus mean roundness for 21 observations 
per operator. 


It is often considered that sphericity and roundness tend to 
vary sympathetically. The means of the 21 independent esti- 
mates of sphericity and roundness for each of the 21 grains are 
plotted in figure 8. No obvious trend appears. 

As a matter of interest, if we plot the mean sphericity and 
roundness for each operator, however (fig. 9), a distinct trend 
is shown. Firstly, the approximate slope of the trends suggests 
that operators tend to react in sympathy to visual estimates of 
roundness and sphericity. In other words, when an operator 
estimates high sphericity, he also tends to estimate high round- 
ness and vice versa. 

Secondly, there is a reversal in the case of five of the seven 
operators which shows that on the first and third occasions 
(days 1 and 3) three of the operators measured high sphericity 
and roundness, whereas on the second day these three estimated 
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relatively lower roundness and sphericity on the same grains. 
Two other operators showed the reverse tendency. This is 
presumably connected with an arrangement of the technique; 
some of the operators remarked that there was a tendency to 
compare the first grain with the comparison chart and there- 
after to compare each grain with the preceding one. To counter- 
act this effect it was suggested (by M.A.R.) that the operators 
run through the grains from numbers 1 to 21 the first day, 
reverse the order of grains the second, and return to the 
original order the third. The reversal of the estimates is doubt- 
less connected with this reversal in arrangement. The incon- 
sistency in the reaction of the three operators contrasted with 
that of the two operators reduced the effect of differences due 
to days—or more correctly, due to order of the grains. 

As was noted during the analysis, the differences between 
day means in sphericity were significant at the 5% level (table 
4), but the differences between day means for roundness were 
not significant. Similarly, variation of grains from day to day 
as estimated by the interaction term grains times days was not 
significant; nevertheless, plotting the day means (inset, fig. 9) 
shows again the sympathetic variation of sphericity with 
roundness. 

The moral is obvious, a more extended test of day to day or 
order of treatment of grains is necessary before these items 
can be considered relatively unimportant. In practice it would 
be advisable to randomize the order of treatment of grains on 
each occasion. Furthermore, the tendency among operators to 
record high sphericity along with high roundness and vice versa 
when using visual comparison techniques suggests some care in 
accepting sympathetic variation in sphericity and roundness 
as an attribute of grains. 

An interesting item in the estimation of roundness is empha- 
sized by the relationship between the mean and variance of 
the 21 grains (fig. 4). The regression indicates that the vari- 
ance increases with decreasing mean roundness, which in turn 
implies that operators can estimate with less ambiguity dif- 
ferences in roundness in well-rounded grains and are increasing- 
ly inconsistent as the grain becomes more angular. In other 
words smaller differences can be detected in grains of high 
roundness than in grains of low roundness. This is in complete 
contradiction to the conclusion stated by Pettijohn (1949, p. 
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51), “ ... because the eye can readily distinguish slight dif- 
ferences in roundness when the roundness values are low, but 
cannot make such distinctions when the values are large.” If 
future experiments confirm our findings, then the roundness 
scale suggested by Pettijohn should be reversed and more 
classes included in the higher range of roundness than the lower. 
Theoretical analysis of these problems is always insidiously 
deceptive, but it does appear reasonable to conclude that small 
departures from a smooth rounded surface are easier to detect 
than on a serrated one! It seems abundantly clear to the writers 
that subjective classification, particularly by visual estimates, 
even with comparison scales such as the charts used in this 
experiment, are no guarantee of valid conclusions from the 
observations unless the experiment is suitably designed; it is 
essential to evaluate the different sources of variation and 
establish that the differences observed are associated with real 
differences in the variable estimated and not introduced by 
lack of control of the psychological bias present in every 
observer. 
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FURTHER DATA ON THE SYSTEM 
Na,O—CaO—Si0, 


E. R. SEGNIT 


ABSTRACT. Primary phase fields, liquidus surface, and quintuple points 
have been determined in part of the system Na,O—CaO—SiO, on the 
silica-poor side of the join CaSiO,-Na,SiO,. A new ternary compound, 
probably Na,Ca,Si,O,,, has been described. The solid solution limit of 
Na,O in dicalcium silicate has been shown to approach 3 per cent. The 
significance of these results in relation to Portland cement clinker com- 
ponents has been discussed, and the conclusion reached that most of the 
Na,O in Portland cement is contained in the dicalcium silicate phase. 


INTRODUCTION 


HE high silica portion of the ternary system soda—lime— 

silica is of fundamental importance to the glass techno- 
logist, containing as it does, the area of compositions of common 
glasses. This portion of the system has, therefore, attracted 
much attention. The initial work was that of Morey and Bowen 
(1925), who have thoroughly investigated the phase equilibria 
of the NS—-CS—S' portion of the system. In this they found 
ten primary phase fields, including those of cristobalite, tridy- 
mite, quartz, aCS, BCS, NS, NS,., and those of the three 
ternary compounds NC,S,, NoCS, and NC,S¢, the latter two 
of which had incongruent melting points. A further ternary com- 
pound, NCS, has been described fully, together with the other 
ternary compounds, by Wyckoff and Morey (1926). Thilo 
(1951) mentioned the compound NC,S.. 

The aim of the present work was to extend the knowledge 
of the phase equilibria of this system into the area containing 
less than 50 per cent silica. In the high-lime portion of the 
triangle occur the fields of the cement compounds C,S and 
CS. The side C-CS is part of the well known lime—silica binary 
system (Rankin and Wright, 1915; Greig, 1927) from which 
limits were obtained on this boundary of the fields for lime, 
C,S, C,5,, and aCS, The join NS-CS forms part of the system 
investigated by Morey and Bowen, and gave limits for the 
aCS, NC.S, and N.CS, fields. The more detailed work done 


1 Throughout the paper, formulae will be abbreviated in this manner, 
viz: N= Na,O, C—CaO, S—SiO,. 
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here has led to their boundary between the NCS, and NS fields 
being slightly altered. This, however, has not affected the es- 
sential equilibrium relationships. 


EXPERIMENTAL 


For the preparation of the mixtures, sodium and calcium 
carbonates of analytical reagent quality were used. The source 
of silica was quartz crushed to pass through a 200 mesh sieve, 
washed with hydrochloric acid and then with distilled water. 
A typical analysis of this product was SiO, 99.86%, R,O, 
0.04%, CaO 0.03%, HO 0.05%. The three constituents were 
weighed in the required proportions, fused in a platinum cruci- 
ble, and, when possible, quenched to a glass. Each mixture 
was crushed and fused four times to ensure homogeneity. All 
glasses in this part of the system crystallized and reached 
equilibrium quite readily below their liquidus temperatures. The 
glassy charges could therefore be used directly for the quench- 
ing experiments without the need for preliminary treatment to 
render them crystalline. 

The primary phase fields, boundary lines, quintuple points, 
and liquidus temperatures were determined by the quenching 
method, Small charges, wrapped in platinum foil of thickness 
0.0005 inch, were rapidly quenched in water after holding at 
a constant temperature until equilibrium was reached. The 
resultant material, commonly a glass plus one crystalline phase, 
was examined with the petrological microscope. 

The electric furnaces were wound with platinum—20% 
rhodium alloy wire, and the temperatures regulated to within 
1°C. by a bridge type electronic controller; the furnace winding 
formed one arm of the bridge. The power into the furnace was 
controlled by a continuously working reversible choke instead 
of the more usual method of shunting the heating current 
through an extra resistance. Temperatures were measured with 
platinum-platinum 13% rhodium thermocouples checked against 
calibrated standards, 

Difficulty was experienced through loss of soda by volatiliza 
tion from compositions of high melting points rich in this con- 
stituent. In the earlier part of the work mixtures containing 
over 10 per cent soda were analyzed after preparation, so that 
the starting compositions were accurately known. In the later 
stages mixtures were prepared at temperatures between 1100°C. 
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and 1150°C., when the vapor pressure of the Na,O was too 
low to cause any appreciable loss. The losses during quenching 
experiments were not easily determined. Many quenches, how- 
ever, were carried out in small buckets containing sufficient 
material for a soda determination. The only quenches seriously 
affected were in those parts of the NCS, NC,S, and C,S fields 
outside the 1300°C. isotherm and containing over 15 per cent 
Na.O. For unanalyzed quenches in this region, allowance was 
made for the final composition by assuming losses similar to 
those which were determined. Throughout the investigation, 
heating times were kept as short as equilibrium requirements 
would allow. 


RESULTS 


In the part of the system investigated nine primary phase 
fields were encountered. These were of the compounds pseudowol- 
lastonite (aCS), dicalcium silicate (C,S), C,S,, NS, NC.S,, 
N.CS,, NCS, NC,S, and N,C,S,. As the first seven of these com- 
pounds have been fully described in the literature, only observa- 
tions pertinent to the present work will be recorded for them. 

C,S: In the higher temperature quenches, charges poor in 
soda gave rise to small individuals, almost spherical in shape, 
which became oval or slightly elongated in melts of high soda 
concentration. In all cases the crystals appeared to be the 
stabilized a-form. 

aCS: Large, well-formed individuals of high birefringence. 
Hexagonal outlines were frequently seen. 

C,S,: Fairly well-formed prismatic crystals of moderate re- 
fractive index and low birefringence. 

NC,S,: Well-formed pseudo-cubic crystals of square outline 
and low birefringence. 

N,.CS,: Small isotropic octahedra. 

NCS: Oval-shaped isotropic crystals. 

NS: Needles or lath-like crystals of moderate birefringence. 

NC,S,: This compound has not previously been described, 
although it is mentioned by Thilo (1951) in connection with 
hydrated soda-lime silicates. In these studies it appeared char- 
acteristically as lath-shaped crystals of low birefringence. Its 
refractive indices are close to 1.665, birefringence 0.003, and 
optical character biaxial with high optic axial angle. The 
extinction of the lathes is inclined, and the crystals are probably 


Further Data on the System Na,O—CaO—Si0, 


SiO, 


80 No,0 40 30 
20 


Fig. 1. The system Na,O—CaO—SiO,. Primary phase fields, composition 
triangles, and compositions studied. (The area NS-CS-SiO, after Morey 
and Bowen, 1925.) 


monoclinic. Many crystals were twinned on a multiple law in 
a similar manner to plagioclase. 

N,C,S,: This previously undescribed compound was first 
seen in melts of high soda concentration as small rectangular 
crystals of very low birefringence and high refractive index. 
Near the liquidus temperatures of melts in which it occurred 
as a primary phase it formed large idiomorphic flat crystals 
showing what appeared to be two large basal planes with two 
sets of dome faces joining them. Its refractive indices are very 
close to 1.620, birefringence 0.002, biaxial +-ve with high optic 
axial angle. The spacings (in A) of lines observed in powder 
photographs of the last two compounds are given in appendix 2. 
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Fig. 2. The system Na,O—CaQ—SiO,. Primary phase fields and liquidus 
isotherms. (The area NS-CS-SiO, after Morey and Bowen, 1925.) 


The results of the more important quenching experiments are 
set out in appendix 1. The phase diagrams showing liquidus 
isotherms, composition triangles, phase boundaries and com- 
positions studied are illustrated in figures 1 and 2. These dia- 
grams extend the work of Morey and Bowen, which is embraced 
by the triangle NS-CS-S, into areas richer in soda and lime. 


DISCUSSION 


Dicalcium silicate.—Dicalcium silicate (Ca,SiO,) has at- 


tracted a great deal of attention due to its presence in a variety 


of industrial materials (e.g., cements, slags, refractories), and 
because of its complex solid phase changes (Bredig, 1943, 
1950). The importance of this compound in the present system 
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lies mainly in its capacity to form solid solutions with small 
amounts of soda, and hence its relationship to the dicalcium 
silicate in high-alkali Portland cement clinkers. This solid solu- 
tion has been recorded earlier by Greene (1944) and Newman 
and Wells (1946). The former prepared samples containing 
Fe,0, or Al,O, in addition to soda, and a dicalcium silicate 
containing 4.2 per cent Na,O and 4.6 per cent Fe,O, was 
identified. The latter authors state that “the limit of solubility 
of Na,O in a-2 CaO°SiO, is uncertain, but probably exceeds 
1.3 per cent.” 

By measuring refractive indices of rapidly cooled prepara- 
tions close to C,S in the soda—lime-silica system, further data 
have now been obtained on the extent of these solid solutions. 
The refractive indices and birefringence of these samples are 
recorded in table 1. It was thought that the C,S was in the a 
(hexagonal) form. The samples were prepared at 1400°C. and 
cooled rapidly. This temperature was probably above the 
a—a’ inversion temperature for the solid solution, and neither 
exsolution effects nor twinning which might indicate inversion 
were seen. 

The measurements in table 1 show that the birefringence 
and refractive indices decrease as the soda content increases at 
least up to 2.7 per cent. The maximum possible soda content 


TABLE 1 


Refractive Indices and Birefringence of Dicalcium 
Silicate Solid Solutions 


Composition -—— Refractive —— 
Na,O CaO SiO. Index Birefringence 


65.1 34.9 1.716 
62.0 38.0 1.703 
64.0 36.0 1.698 
66.0 34.0 1.712 
62.4 37.1 1.712 
63.2 36.0 1.709 
1. 
1 


0.20 
.028 
.029 
021 
020 
021 
015 
O19 
O19 
O15 
O15 


66.4 32.1 708 
63.3 35.0 705 
63.6 34.3 712 
62.0 35.9 
65.5 32.3 1.710 
2.7 62.3 35.0 1.698 


*dicalcium silicate composition. 


736 
731 
134 
733 
7129 
729 
.720 
731 
.728 
.725 
713 
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is probably in the region of 3 per cent, and from the limited 
data available, seems to be greatest in the directions of the 
compounds NCS and NC.S,. The change in birefringence was 
well marked in the quenching experiments, where a distinct 
lowering of the interference colors of crystals was seen in 
primary phase C,S quenched from melts containing quite small 
quantities of soda. 

It seemed probable that at least some, and possibly a con- 
siderable part of the soda in Portland cement occurs in the 
dicalcium silicate. Ordinary Portland cement clinker generally 
contains about 25 per cent C,S, so that ideally some 0.75 
per cent Na,O could be taken up in this compound. Any soda 
above this upper limit would have to be accounted for in other 
phases, such as the glassy interstitial material. The dicalcium 
silicate in three laboratory-prepared Portland cement clinkers 
of varying soda content was examined optically. The refractive 
indices are recorded in table 2, from which it is evident that 
the refractive index decreases with the soda content of the 
clinker. The only other constituents of the cement which might 
enter into the C,S structure and affect its properties are Al,O, 
and Fe,O,. Of these, the former would have little effect (un- 
published work by the author on other materials containing 
Al,O,), and the latter would tend to increase the refractive 
index. The decrease in refractive index shown in table 2 can 
therefore be assigned with reasonable certainty to the effect 
of soda. 

The optical properties tabulated in table 1 also show that 
the refractive indices are decreased and birefringence increased 
in soda-free compositions which contain an excess of silica over 
the stoichiometric composition. This may be due to a small 


TABLE 2 


Refractive Indices of Dicaleium Silicate from Laboratory 
Portland Cement Clinker 


Clinker Na,O Content -— Refractive Indices — 
Number of Clinker of C,S from Clinker 


26 B 2.4 er 1.71 
27 B 1.0 oer 1.72 
22 B < 9.1 1.715 1.73 
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amount of solid solution with aCS, which has a birefringence of 
.044 compared with .020 for C,S. Newman and Wells showed 
that the solid solution limit of SiO, in dicalcium silicate is 
0.5 per cent (approximately 1 per cent CS). It may be noted 
that excess CaO or SiO, in the C,S inhibited to a large extent 
the inversion to the y form. 

Equilibrium relationships.—These are summarized in figures 
1 and 2 which show the primary phase fields, liquidus isotherms 
and composition triangles. Thirteen new boundary lines have 
been located in part or in full, and nine quintuple points. Of 
the latter, four were eutectics. Parts of nine primary phase 
fields have been mapped. Crystallization paths and other data 
may be interpreted from figures 1 and 2. A few points, however, 
deserve further mention. 

The melting point of the compound NC,S, (Na,Ca,Si,0;) 
was not determined with great accuracy owing to the high 
vapor pressure of soda at these temperatures. A sample of the 
compound quenched from 1446°C. was almost entirely glass, 
so that, allowing for a small loss of soda during heating, the 
melting point of the pure compound would be a little higher 
than this, and close to 1450°C. Its primary phase field is 
bordered by those of aCS, NC,S,;, C,5, C,5. and NCS and 
four of the invariant points (B, C, D and E) at the corners 
of its field are eutectics. 

The compound N,C,S, (Na,Ca,Si,0,;) melts incongruently 
at 1125°C. It has a small primary phase field FGHI some dis- 
tance from its composition point, which was located only with 
some difficulty. As there were no maxima on any of the bound- 
aries of its field, it therefore probably lay to the right of the 
lines NS-F and F-NCS. Several possible compounds in this area 
were prepared by repeated sintering at 1000°C. and grinding, 
and it was found that the above composition satisfied the re- 
quirements of the compound. The product was composed almost 
entirely of the crystals obtained as a primary phase in the 
area GHI together with a small amount of glass.of refractive 
index 1.55. The only simple molecular ratios of the three oxides 
richer in lime in the possible area are N,;C,S, and N.C,S;. 
On preparation the former gave a mixture of NCS and the 
new compound. The latter, while containing a large proportion 
of the new compound, gave certain additional strong lines in 
a powder photograph, which did not appear in the pattern of 
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the 4:3: 5 compound. The invariant point F is a little unusual 
in that the temperature falls along both boundaries FG and FI, 
and is therefore a “falling reaction point.” No eutectic was 
found in the soda-rich part of the area investigated. 

The presence of a C,S, field is indicated by the appearance 
of a small liquidus for C,S, in the lime-silica binary system. 
A series of quenches on charges of the composition 5% Na,O, 
50% CaO, 45% SiO, which lies near the edge of the C,S field, 
shows that a long, narrow field extends to nearly 9% Na,O. 
At the reaction point A C,S, is in equilibrium with C,S, 
NC.S, and liquid, and at the eutectic B is in equilibrium with 
aCS, NC,.S, and liquid. 

In table 3 are given the composition and temperatures of the 
nine invariant points. 


TaBLe 8 


Ternary Invariant Points 


——Composition— 
Na,O CaO 
Point Phases Type % 


C,S-C,S.-NC,S, Reaction . AL 1280 
a CS-C,S,-NC,S, Eutectic 1270 
NC,S,-a CS-NC,S, Eutectic 9.! 1255 
NC,S,-NCS-NC,S, Eutectic 4. 82. 1260 
C,S-NCS-NC,S, Eutectic 3. 37. 1440 
NC.S,-NCS-N,C,S,; Reaction A 1130 
NC,S,-N.CS,-N,C,S; Reaction 40.0 1120 
NS-N.CS,-N,C,S, Reaction 45.0 990 
NS-NCS-N,C,S, Reaction 48.0 950 
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APPENDIX 1 


Results of Quenching Experiments 


--Composition Wt. Temperature Phases 
Na,O CaO SiO, *C. Observed 


13.0 39.0 48.0 1286 Glass 
1281 Glass, rare NC,S, 
field 1275 Glass, NC,S, 


1287 Glass 
1280 Glass, NC,S, 
1255 NC,S,, NC,S,, tr. glass. 


1270 Glass, NC,S, 
1260 NC,S,, NC,S,, glass 


1280 Glass 
1266 Glass, NC,S, 


1285 Glass 
1283 NC,S,, rare glass 


1274 Glass 
1269 Glass, NC,5, 


1288 Glass, NC,S, 


1240 Glass 
1236 Glass, NC,S, 


~ 
20.0 30.0 50.0 
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Results of Quenching Experiments 


-—Composition Wt.%— Temperature Phases 
Na,O CaO SiO, A Observed 


23.0 80.4 46.6 Glass 
Glass, NC,S, 
NC,S,, glass 
NC,S,, NCS, tr. glass 


Glass 
Glass, NC,S, 


Glass, rare NC,S, 
Glass, NC,S,, NCS. 


Glass 
Glass, NC,S, 


Glass 
Glass, tr. NC,S, 


Glass 
Glass, NC,S, 


Glass 
Glass, NC,S, 


Glass 
Glass, N,CS, 


Glass, rare N,CS, 
Glass, N,CS,, NS 
NS, N,C;S,, N,CS, 


The 24. 35. 40. 4! Glass 
NCS Glass, NCS 
field 
Glass 
Glass, NCS 


Glass 
Glass, NCS 


Glass 
Glass, NCS 


Glass 
Glass, rare NCS 


The 
field 
23.0 82.0 45.0 1266 P| 
1258 
24.4 30.2 45.4 1269 Po 
1258 
29.6 25.1 45.3 1209 
1199 
81.0 26.7 42.3 1234 P| 
1229 
34.4 20.2 45.4 1149 po 
1141 
38.3 18.2 43.5 1115 P| 
1101 
The 38.2 15.4 45.4 1128 P| 
N,CS, 1117 
field 
44.5 10.0 45.5 1010 
988 
961 
29.5 31.2 89.3 1348 P| 
1340 
29.7 30.1 40.2 1379 P| 
1370 
81.0 31.3 87.7 1452 P| 
1429 
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--Composition Wt.%— Temperature Phases 
Na,O CaO SiO, << Observed 


32.8 27.0 40.2 Glass 
Glass, NCS 


33.7 25.5 0. Glass 
Glass, rare NCS 


23.5 ‘ Glass, rare NCS 


20.4 : ‘ Glass, rare NCS & NC,S, 
NCS, NC,S, 


Glass 
Glass, rare NCS 


Glass, rare NCS 
Glass, NCS 


Glass, NCS 
Glass, NCS, N,C,S, 


Glass 
Glass, NCS 
Glass, NCS, N,C,S; 


Glass, rare NCS 
Glass, NCS, N,C,S, 


Glass, NCS 
Glass, NCS, N,C,S, 


Glass 

Glass, NCS 

Glass, NCS 

NCS, NS, N,C,S,, glass 


Glass, rare NCS 
Glass, NCS 


Glass, rare NCS 
Glass, NCS, rare NS 
Glass, NCS, NS 


Glass, NC,S, 
Glass, NC,S,, rare CS 


597 
NCS 
field 
4.60 
89.1 18.3 42.6 1167 Po 
1157 
39.8 20.1 40.1 1317 P| 
1290 
40.0 17.5 42.5 1145 
1124 
41.3 18.7 40.0 1159 
1139 
1096 
42.2 17.0 40.8 1120 Po 
1091 
42.9 15.5 41.6 1065 Po 
1005 
44.5 15.1 40.4 1136 
1125 
1003 
949 
47.0 13.0 40.0 1100 P| 
1073 
49.0 10.0 41.0 942 
929 
921 
The 11.0 43.0 46.0 1289 P| 
NC,S; 1268 
field 
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Results of Quenching Experiments 
Composition Wt. of, Temperature Phases 
Na.O CaO SiO, A Observed 
11.0 44.0 45.0 Glass 

Glass, NC.S,,C.S 


Glass 
Glass, NC,S,,C,5S; 


Glass 
NC,S,, NC,S,, glass 


Glass 
Glass, NC,S, 


Glass 
Glass, NC.S, 


1303 Glass 
1300 Glass, rare NC, 


1384 Glass, rare NC.S 


1389 Glass, rare NC,S 
1366 Glass, NC,S, 


1315 Glass 
1307 Glass, NC,S, 


1446 Glass, rare NC,S, 
1442 Glass, NC.S, 
1435 NC,S, 


1370 Glass 

1353 Glass, NC,S, 
Glass 
Glass, NC,S., 
Glass, NCS 


Glass, rare NC,S, 


The 43.1 : Glass, rare N,C,S, & 
N,C,S, NC,S, 
field 


1092 Glass, rare N,C,S, 


1092 Glass, rare N,C,S, 


we 
The 
NC,S, 
field 
1267 
15.0 38.0 47.0 1268 
1261 
15.0 40.0 45.0 1347 P| 
1332 
17.0 37.0 46.0 1310 P| 
1301 
S, 
18.0 38.0 44.0 
22.0 35.0 43.0 
23.0 36.0 41.0 
24.0 33.0 43.0 1327 P| 
410 160 43.0 
41.5 15.0 43.5 
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Results of Quenching Experiments 


Composition Wt.%— Temperature Phases 
Na,O CaO SiO, “C. Observed 


43.2 15.0 41.8 1046 Glass 
1038 Glass, N,C,S, 


1050 Glass, rare N,C,S, 
1045 Glass, N,C,S, 
Glass, N,C;S,, NS 


Glass, rare N,C,S, 
Glass, N,C,S, 
Glass, N,C;S,, NS 


Glass 
Glass, rare N,C,S, & NS 
Glass, N,C,S,, NS 


Glass 
Glass, N,C,S, 
Glass, N,C,S,, NS, NCS 


Glass 
Glass, NS 


Glass 
Glass, NS 


Glass 
Glass, rare NS 
Glass, NS, N,C,S, 


Glass 
Glass, NS 
Glass, NS, N,C,S, 


Glass 
Glass, NS, NCS 


1065 Glass 
1058 Glass, NS, rare NCS 


1069 Glass 

1061 Glass, NS 

1438 Glass 

1426 Glass, C,S 

1400 Glass, C,S,,a CS 


The 
field 
44.0 13.3 42.7 1024 
1014 
974 
45.6 10.3 44.1 989 
981 
968 
47.0 12.0 41.0 1035 
1005 
950 
The 46.0 10.0 44.0 983 
NS 979 
field 
46.5 8.4 45.1 1012 
1006 
47.0 10.0 43.0 976 
971 
959 
47.2 7.1 45.7 1031 
1019 
984 
48.0 10.0 42.0 967 
963 
| 
51.0 2.0 47.0 
The 4.0 51.0 45.0 
field 
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Results of Quenching Experiments 


-—-Composition Wt.%—, ‘Temperature Phases 
Na,O CaO SiO, << Observed 
The 4.9 55.1 40.0 Glass, C,S 
field 5.0 45.0 Glass, rare C,S 
Glass, C,S 
Glass, C,S,C,S, 
Glass, C,S,,CS 
C,S., a CS, glass 


Glass, C,S 


Glass, rare C,S 
Glass, C,S, NC,S, 


Glass, rare C.S 
Glass, C,S 


Glass 
Glass, C,S 


Glass, C,S 
C,S, glass 


C,S, NCS, glass 


C,S, rare NCS, glass 


Glass 
Glass, a CS 


Glass, a CS 


Glass,a CS 


7.9 51.2 40.9 1509 FY 
9.0 45.0 46.0 1324 P| 
1307 
10.0 45.0 45.0 1341 P| 
1338 
17.5 42.3 40.2 1453 P| 
1449 
20.0 40.0 40.0 1460 
1442 
23.2 46.1 80.7 1580 PF 
27.1 428 30.1 1585 Po 
The 9.0 44.0 47.0 1340 po 
1339 
field 
12.0 37.0 51.0 1326 Fe 
12.1 388.7 492 1296 
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APPENDIX 2 
X-ray Diffraction Lines of NC2S» and N4C3S5 
din A 


—- NC,S, ——————_, 
d 
1.836 
1.814 
1.717 
1.619 
1.599 
1.586 
1.554 


< 
= 


WwW 

WwW 
WwW 
WwW 
WwW 
WwW 
WwW 
WwW 
WwW 


= 
< 
= 


< 
< 
=< 


< 
= 
= 


d I d I d I 
4.04 W 6.71 Vw 2.171 VW 
3.86 W 5.20 VW 2.030 § 
3.57 M 4.76 Vw 1.972 Ww 
3.29 WwW 4.41 Vw 1.928 Vw 
8.15 Ww 
2.90 W- 
2.76 vs 
2.71 vs 1.536 
2.61 s 1.514 
2.59 s 1.371 
2.50 Vw 1.350 
2.41 WwW 1.327 
2.31 WwW 1.300 
2.27 WwW 1.289 
2.22 WwW 1.276 
2.18 Ww 1.263 
2.11 Vw 1.216 
2.04 § 1.116 
1.966 VW 1.085 
1.922 VW 1.018 
1.866 W 
2.251 
2.219 
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RUIN-MARBLE: 
A STYLOLITIC FORMATION 


B. M. SHAUB 


ABSTRACT. Ruin-marble, a banded calcareous material, is described as 
a structure related to stylolites in which the many column-like projections 
are cross sections of more or less small tabular fault blocks which have 
originated through contraction and pressure phenomena accompanying 
the dewatering of the calcareous sediments. This structure indicates the 
possibility that stylolites can possibly develop under certain conditions 
without the presence of a clay seam. There is no evidence that solution 
has occurred at any time during the development of the ruin-marble. The 
more highly ferruginous layers remained more plastic than the calcareous 
layers and were occasionally pressed into the fractures and other places 
where there was a sufficient release of pressure. 


eaten is known chiefly from the environs of 
Florence, Italy, and also from Klosterneuburg, Austria. 
It is usually displayed in mineral collections along with such 
materials as landscape agate and other chalcedonic structures 
which contain included dendritic configurations of various kinds. 
In their manner of origin and their appearances there is, how- 
ever, a very great difference. Little information appears to be 
available in the literature regarding ruin-marble. Dana (1889, 
p. 679) states that “ruin-marbie is a kind of compact calcareous 
marl, showing when polished, pictures of fortifications, temples, 
etc. in ruins due to infiltration of oxide of iron;...” Merrill 
(1910, p. 340) describes ruin-marble as a “very compact yellow- 
ish or drab limestone, the beds of which appear to have been 
fractured in every conceivable direction by geological agencies, 
after which the resultant fragments have become cemented by a 
calcareous or ferruginous cement. The rock is therefore really 
a breccia, although the proportional amount of cement is very 
small, and the actual displacement of the various particles but 
slight. When cut and polished the slabs have somewhat the 
appearance of mosaics representing the ruins of ancient castles 
or other structures. Hence the name of ‘ruin-marble’.” 

Because of the lack of suitable material to conduct a more 
comprehensive study of the structure, a tentative description 
will be given in order to bring again this unusual structure 
to the attention of those interested in sedimentary processes and 
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the structures developed therein during the period of dewatering 
and consolidation. 

Examination of the photographs (plate 1) shows the pres- 
ence of certain prominent dark layers, rich in iron oxide, 
which carry across a number of the “columns” in the same 
stratigraphic sequence, while in others the same prominent dark 
bands appear to be missing and a new sequence of banding ap- 
pears. One of the prominent features is the pronounced wedg- 
ing of the columns, showing not only that the relative move- 
ments have been up and down with lateral displacements to the 
right or left but also there have been similar displacements to 
the front and rear of the sections shown in the photographs. 
A small piece of material from Klosterneuburg, Austria, was 
polished on three planes at approximately right angles, and in 
each plane components of the displacements were clearly shown 
to have occurred in two directions in each plane. The bound- 
aries between the various parts are in reality minute fault 
planes, the continuation of which can often be traced into the 
“overlying” unbanded material, according to the orientation 
of the illustrations. In the few pieces examined there were no 
open cavities along these fault planes. This is probably due 
to the retention by the material of a sufficiently high degree 
of plasticity to prohibit the formation of cavities at this stage; 
or any cavities encountered during the preparation of the mate- 
rial for market would probably have been discarded. 

It is quite evident from the illustrations that the iron oxide 
bands retained a higher degree of plasticity during the period 
of the pronounced adjustment, for they were frequently squeezed 
into bulges and along the many small fault planes as shown 
at a, b, c, d, and e in figure 1 and at m, n, o, and p in figure 2. 
There are also instances of slight bending of the bands along 
the contacts, but in no more than a few instances is the bending 
prominent. This shows that there was no heavy lateral pressure 
along the small fault planes at the time of displacement and that 
the development progressed slowly under very low differential 
pressures continually regenerated here and there through the 
structure during the dewatering process. 

It is, indeed, certain that there has not been any solution of 
the calcareous material before, during, or after the formation 
of the structure, for there is not the slightest indication that 
there has been any removal of material by solution. 
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The origin of the structure appears to be a simple straight- 
forward process which consisted of shifting small masses of the 
banded calcareous material from place to place through the 
combined action of the contraction and slight differential pres- 
sures during the period of contraction during dewatering. The 
pressures operating were a combination of the negative contrac- 
tion forces and the static pressure due to the overlying material. 
The contractions produced a lateral relief of pressure which 
permitted the interfingering of the various columns as they 
now appear. They are in reality quite irregular tabular fault 
blocks of minute proportions. While the process is a simple 
one, the relative movements of the pieces are very complex and 
intricately interrelated to one another. Sections may appear to 
be missing but they have been displaced to the front or rear 
of the section in view. 

This structure appears to be closely related to stylolitic 
development. Its unique characteristics have been inherited 
through its peculiar sedimentary environment. There are, how- 
ever, some striking differences between ruin-marble and the 
various stylolitic seams and columns. One feature of the ruin- 
marble which first attracts attention is the clearcut sharp 
angular outlines of the blocks or sections which are in marked 
contrast to the more or less continuous nature of the material 
on the two sides of a normal stylolite seam. It is the sharp 
and pronounced banding which shows the relative movements 
of the blocks in the lower part of the illustration. However, 
the non-banded part above has probably undergone a similar 
adjustment as the fracture planes continue into this material. 
A few waterlines appear to cross all fault blocks at the same 
level. A similar adjustment, but probably of a more uniform, 
continuous nature occurs on the two sides of most of the normal 
stylolite seams, for the Tennessee marble carries a very great 
number of indications of plastic adjustment in the mass of the 
marble on the two sides of the seams. 

The striking difference between the ruin-marble and the 
normal stylolite seam is the absence of a clay seam which usually 
accompanies and is an essential to the normal stylolite seam, 
although the clay partings may be almost infinitesimally thin 
at times. The absence of clay seams here shows that under 
certain conditions they are probably non-essential in the devel- 
opment of stylolites. The movement of the residual water into 


Ruin-marble from Florence, Italy. 


Ruin-marble showing a very intricate rearrangement of the 
banded lavers of calcareous material. Locality unknown. 
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and along a stylolite seam after its development by plastic 
flow can deposit the minute amount of colloidal clay-like material 
often found along such seams. 

Acknowledgment.—The writer wishes to express his apprecia- 
tion to Mr. Howard E. Norris, curator of the Henry Miller 
mineral collection on exhibit at the Loomis School at East 
Windsor, Connecticut, for the privilege of examining and photo- 
graphing the ruin-marble specimen from Florence, Italy. 
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DOWNSTREAM CHANGE OF VELOCITY IN 
RIVERS 


LUNA B. LEOPOLD 


ABSTRACT. Because river slope generally decreases in a downstream 
direction, it is generally supposed that velocity of flow also decreases 
downstream. Analysis of some of the large number of velocity measure- 
ments made at stream-gaging stations demonstrates that mean velocity 
generally tends to increase downstream. Although there are many reaches 
in nearly all rivers where mean velocity decreases downstream, the general 
tendency for conservation or for downstream increase was found in all 
data studied. 


Computations of bed velocity indicate that this parameter also tends 
to increase downstream. 

Near the streambed, shear in the vertical profile of velocity (rate of 
decrease of velocity with depth) tends to decrease downstream. This down- 
valley decrease of shear implies decreasing competence downstream. 


GENERAL STATEMENT 


ECAUSE rivers tend to be steep in their headwaters and 

flow at flat gradients near their mouths, it appears logical 

to suppose that the stream velocities are great in the headwaters 
and decrease downstream. 

Yet, the engineer who actually measures stream velocities for 
purposes of computing the water discharge uses the same 
current meter for a measurement in a small brook as in the 
Mississippi River. The velocities that he measures are not suf- 
ficiently different in small headwaters and downstream reaches 
to furnish any obvious indication of a definite pattern of 
organization in the data. 

The purpose of the present paper is to organize some of the 
large amount of available data on river velocities so that the 
downstream change may be evaluated. 


DATA ON RIVER VELOCITY 


There are more than 6000 gaging stations in the United 


States. Records on many streams are more than 50 years in 
length. At each station periodic measurements of discharge are 
made by use of a current meter. The discharge on other days 
is determined indirectly from the recording chart or daily read- 
ings of stage (water surface elevation). 
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When a river measurement is made by current meter, the 
water speed is measured at about 20 to 30 verticals across 
the section and usually at two depths in each vertical. From 
these direct measurements the form of the river cross section 
may be plotted, as well as an approximate pattern of the 
vertical and horizontal distribution of velocity. 

In stream-gaging practice, the mean velocity is not a simple 
average of the velocity measurements made at various positions 
across the stream. The current meter placed in a given position 
in the cross section measures the downstream movement of the 
water in a given interval of time. The metering thereby defines 
a prism of water whose width and depth are known and whose 
downstream dimension is represented by the measured move- 
ment. The volume of the prism that passed the meter during 
the measured time interval represents the discharge in units 
of volume per unit time as denoted by the identity, 


Q=aV 


where Q is discharge in cubic feet per second, a is cross-sectional 
area in square feet, and V is velocity in feet per second. The 
total discharge past the cross section is the sum of the dis- 
charges represented by the individual prisms. 

The total discharge having been obtained in this way, the 
mean velocity used in hydraulic practice is an expression of the 
velocity of a single rectangular prism having the same cross- 
sectional area as the stream. This mean velocity is defined as 
the quotient of the discharge in the cross section divided by 
the cross-sectional area of flowing water. Arithmetically the 
mean velocity so defined is equivalent to a weighted mean in 
which the velocities of the individual prisms are weighted by the 
cross-sectional area of the respective prisms. 

Owing to the distribution of the velocity with depth, the mean 
velocity is closely equal to that existing at a position 0.6 of 
the distance from the water surface to the streambed. 

It is recognized (Rubey, 1938, p. 129) that the mean 
velocity in the cross section is not the most meaningful velocity 
with respect to sediment transport. Direct measurements of the 
bed velocity, however it may be defined, are not possible at 
present. On the other hand, a large fund of measurements of 
mean velocity exists. It appears instructive to examine the 
downstream change of that velocity parameter for which meas- 
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urements exist, recognizing that the parameter, mean velocity, 
is not ideal for discussing geomorphic progresses. In a later 
section, downstream changes of a computed bed velocity will 
be presented. 

As indicated previously, velocity in a stream is a conserva- 
tive quantity. Even the maximum point velocity (maximum at 
any point in the cross section) does not ordinarily greatly 
exceed the mean velocity for the cross section. An unpublished 
study by the Geological Survey of 2950 measurements of maxi- 
mum point velocity values from a variety of rivers showed a 
median value of 4.11 ft per sec., mean of 4.84 ft per sec. and 
less than one per cent of the total exceeded 13 ft per sec. 

The largest value of maximum point velocity in a natural 
river channel ever measured by stream-gaging personnel of 
the U. S. Geological Survey was about 22 ft per sec. Approxi- 
mately this value was measured on several occasions: Brazos 
River at Waco, Texas, Sept. 27, 1936; Canadian River near 
Newcastle, Okla., May 4, 1941; New River at Caperton, W. Va., 
Aug. 1940; Potomac River at Washington, D. C., May 14, 1932. 


FLOW DURATION AND FLOOD FREQUENCY 


In order to compare the mean velocity in a small headwater 
tributary with that in a large trunk stream it is necessary 
to establish a basis of comparability. It would be meaningless 
to compare the velocity during a flood flow at one point on 
a river with that of a low flow at some other point. The two 
points to be compared must both be experiencing flood flow 
or must both be experiencing low flow. The relative magnitude 
of flow is discussed in hydrologic practice by the use of the 
flood-frequency concept. If the mean flow for each day of record 
at a station is tallied in categories of size, the number of days 
on which each category of size is equaled or exceeded may be 
counted. The number of days in each category may be ac- 
cumulated from the largest category to the smallest, and divided 
by the total number of days. The resulting quotients represent 
the percentage of the total days each given discharge is equaled 
or exceeded. When the discharge equaled or exceeded is plotted 
against per cent of time the resulting curve, called the “flow- 
duration curve,” is directly analogous to the plot used in 
sedimentation work to describe the size-frequency distribution 
of a sediment sample. 
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The velocity at two stations along a river system may be 
logically compared when the discharge at each represents the 
same percentage of time equaled or exceeded. 

The flow-duration curve adequately describes the per cent 
of time a given mean daily flow is exceeded. A flood peak, 
however, may last only a few hours and therefore would not 
be represented by the mean flow for the day on which it occurs. 
The frequency of flood peaks, therefore, is usually determined 
by the use of a series of data representing the momentary 
peak-discharge values for each year of record. This series is 
arranged in order of magnitude, and plotted in terms of the 
recurrence interval in years versus the peak discharge. 

There are minor variations in methods of preparing flow- 
duration and flood-frequency curves, and procedures for anal- 
ysis are adequately described in the literature (Jarvis, 1936, 
pp. 68-85; Linsley, Kohler, and Paulhus, 1949, pp. 182-198), 
but the general principles underlying all methods are essen- 
tially similar. In the analysis contained herein, the procedures 
for obtaining flow-duration and flood-frequency curves were 
those in current use in the U. S. Geological Survey. 


MEAN VELOCITY AT A GIVEN RIVER CROSS SECTION 


The mean velocity in a river cross section tends to increase 
with discharge. This becomes evident from a consideration of 
a relation widely used in hydraulic practice, the Manning equa- 
tion, expressing the balance of frictional losses (which are 
dependent on velocity, hydraulic radius, and roughness) and the 
slope of the energy grade line. 


V=15R% 


where V = mean velocity 


area of cross section 


R= = hydraulic radius 


wetted perimeter 


s = slope of the energy grade line, equal to slope of water 
surface for steady, uniform flow 


n = roughness parameter 


For most natural river channels the width is great compared 
with depth and the hydraulic radius, R, is approximately equal 
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to the mean depth, d. The discharge, Q, is expressed as the pro- 
duct of mean velocity and cross-sectional area, 


With increasing discharge in most natural river cross sec- 
tions, the water-surface slope remains nearly constant and 
roughness tends to decrease. Thus velocity tends to increase 
with increasing depth and therefore with increasing discharge. 

Figure 1 presents a typical relation of mean velocity to dis- 
charge in a river cross section. The data are representative 
of the large volume of observational material collected at stream- 
gaging stations by the use of current meter. Such data are 
available in the files of the Water Resources Division, U. S. 
Geological Survey. For the station pictured, the mean annual 
discharge during a 9-year period of record was 881 cfs, at 
which the corresponding velocity is 3.2 ft per sec. At any given 
discharge a corresponding mean velocity may be read off the 
graph. 


CHANGE OF MEAN VELOCITY DOWNSTREAM 


It has been demonstrated (Leopold and Maddock, 1953) that 
the velocity corresponding to mean discharge tends to increase 
with increasing discharge downstream. 

The reason for the downstream increase of velocity with in- 
creasing discharge has been previously discussed in detail (Leo- 
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Fig. 1. Relation of mean velocity to discharge in a given river cross 
section, Stillwater River near Absarokee, Montana. 
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pold and Maddock, 1953, pp. 15, 42, and 43). The explana- 
tion can be summarized by saying that the downstream increase 
of depth overcompensates for the decrease in slope. 

In that study, the river data considered in analyzing the 
downstream change of velocity did not include enough high 
mountain streams of very steep gradient and in which caliber 
of load rapidly decreased downstream. Moreover, downstream 
changes in bed velocity were not discussed. The present paper 
reviews the evidence previously discussed, presents data for 
streams representing a more stringent test of the generalization, 
and discusses similar relations involving bed velocity. 

The mean annual discharge at a given station is the arithmetic 
mean of all the daily discharge values during the period of 
record. In most river systems, the mean discharge increases 
downstream owing to the increase in drainage area in the down- 
stream direction. As a general rule, mean discharge increases 
directly with some power (near unity) of drainage area in 
a given basin. Though there are many basins, particularly in 
the western United States, in which mean discharge decreases 
owing to irrigation diversion and to channel losses, in most 
basins there is an increase downstream. 

The mean discharge at most points along a river system 
occurs somewhat less often than half the time. The per cent 
of time the mean discharge is equaled or exceeded is not constant 
but is sufficiently conservative to allow a comparison of up- 
stream and downstream points. Because it is a parameter readily 
available for a large number of gaging stations, its use allows 
the comparison of a considerable number of stations. The com- 
putation of complete flow-duration and flood-frequency curves 
has been completed for a smaller number of selected stations 
chosen for a more critical evaluation. 

The mean velocity at a number of stations along a river 
length will be compared at discharges corresponding to their 
respective mean annual values. Figure 2a presents a sample of 
this relation for a small portion of the upper Missouri River 
basin and a few of stations on the Missouri and Mississippi 
rivers. 

Each point in figure 2a represents a particular gaging station 
and its abscissa value corresponds to its mean annual discharge. 
Within each tributary basin, the points are numbered in se- 
quence from downstream toward the headwaters. In nearly all 
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basins the mean annual discharge increases downstream but 
there are examples, such as the Greybull River, in which the 
mean discharge at a downstream station, Basin, Wyoming 
(point 23), is less than at Meeteetsee farther upstream (point 
24), owing to irrigation diversion and channel losses. 

Stations in the Yellowstone and Bighorn River basins chosen 
to represent a headwater area have velocities at mean annual 
discharge of 2 to 4 ft per sec., whereas the lower Missouri 
and Mississippi rivers have corresponding velocities between 
4 and 6 ft per sec. The mean line in the graph of figure 2a 
has been fitted by inspection and is intended merely to call 
attention to the definite tendency of the points to progress 
toward higher velocity values downstream. 


Fig. 2. Change of mean velocity as mean annual discharge increases 
downstream. In 2a above, each point represents a gaging station on Yellow- 
stone River system and downstream. In 2b, mean trend lines for 
various river basins are drawn. 


Identification of Stream-gaging Stations 
in Figure 2A 


- Mississippi R. at Vicksburg, Miss. 
” Memphis, Tenn. 
St. Louis, Mo. 
Hermann, Mo. 
Kansas City, Mo. 
St. Joseph, Mo. 
Pierre, S. D. 
Bismark, N. D. 
Williston, N. D. 
Sidney, Mont. 
Miles City, Mont. 
Billings, Mont. 
. Corwin Springs, Mont. 
. Bighorn R. Custer, Mont. 
‘ St. Xavier, Mont. 
Kane, Wyo. 
Thermopolis, Wyo. 
Manderson, Wyo. 
Crowheart, Wyo. 
Riverton, Wyo. 
Burris, Wyo. 
, Dubois, Wyo. 
. Greybull R. at Basin, Wyo. 
” Meeteetsee, Wyo. 
. Owl Creex nr. Thermopolis, Wyo. 
. Medicine Lodge Cr. nr. Hyattville, Wyo. 
. Gooseberry Cr. nr. Grass Creek, Wyo. 
. N. Fk. Owl Cr. nr. Anchor, Wyo. 


” 


. Missouri R. 
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In any given reach of river the velocity may decrease down- 
stream, as can be seen by a study of figure 2a. Also there is 
a considerable scatter of points in the graph. But for present 
purposes attention is directed to the general trend which is 
toward increasing velocity downstream. 

In figure 2b the trend lines of velocity versus mean annual 
discharge are plotted for a number of rivers chosen to represent 
a variety of physiographic provinces. The general tendency for 
increasing velocity as mean annual discharge increases down- 
stream is well exemplified. 

The mean annual discharge in many streams is equaled or 
exceeded about one out of every 4 days in a long period. 
Inasmuch as the greatest total sediment load is carried by a 
stream during high or flood flows, it is of interest to examine 
the trend of velocity at discharge rates less frequent than the 
mean annual rate. For this purpose figure 3 presents the trend 
lines for discharges equaled or exceeded various percentages of 
time from 50 per cent to 1 per cent, in the Maumee and Scioto 
River basins in Ohio. The two basins were considered together 
in order to be able to utilize a sufficient number of points to 
define the downstream trends. The Maumee and Scioto basins 
were chosen because complete flow-duration curves for stations 
in Ohio are available in the study by Cross and Bernhagen 
(1949). Figure 3 is taken from Leopold and Maddock (1953). 

The scatter of points about the lines represented in figure 
3 was again great, but the general trend of increasing velocity 
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Fig. 3. Relation of mean velocity to daily discharge equaled or exceeded 
various percentages of time, Maumee and Scioto river basins, Ohio. 
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downstream in these basins appears to characterize each of 
the duration percentages analyzed. 

But the Maumee and Scioto basins still do not furnish a 
critical test of how mean velocities vary downstream because 
the topographic relief in those drainage basins is not great. 
Moreover, the materials composing the beds of those rivers 
are relatively fine grained, consisting for the most part of clay, 
silt, and sand. It is logical to suppose that if in any stream 
the velocity tends to decrease downstream, such a stream would 
be one that is very steep in the headwaters and flat downstream. 
Moreover, it would be a stream the bed of which consists of 
coarse gravel in the headwaters, decreasing in grain size to 
sand downstream. Furthermore, the graphs presented earlier 
were for discharges exceeded no less often than 1 per cent of 
the days, or one day in 100. It might be supposed that the more 
extreme floods, which are the most effective in moving coarse 
gravel in the headwaters, would be the flow condition most 
likely to show a tendency for decreasing velocity downstream. 


As an example of these stringent conditions a series of trib- 
utaries of the upper Yellowstone River in Wyoming and 
Montana were selected. The streams chosen include the Grey- 
bull River, Rock Creek, the Clarks Fork, and some of their 
headwater tributaries. These streams drain into the Bighorn 
River, a tributary to the Yellowstone River, or directly into 
the Yellowstone. These headwater tributaries carry coarse gravel 
which decreases in caliber rapidly downstream. The slopes of 
the uppermost tributaries are up to 100 feet per mile. More 
than one tributary to the Yellowstone was selected in order that 
a sufficient number of gaging stations would be included to 
provide an adequate test. 


In accordance with good hydrologic practice, the flood- 
frequency curves for all gaging stations were developed for an 
identical period of record, 1923-1950, inclusive, and the records 
shorter than this were adjusted to the long record. 


The flood data utilized at each station consisted of the peak 
or momentary discharge for each year of available record. 
Flood-frequency curves were developed for recurrence intervals 
of 5, 10, 25, and 50 years. By definition, a flood of 50-year re- 
currence interval means the momentary discharge equaled or 
exceeded on the average once in 50 years over a long period. 
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Fig. 4. Change of mean velocity with increase of discharge downstream 
for flood flows, Yellowstone River Basin and downstream. Upper figure 
shows how river slope decreases downstream. Lower figures show down- 
stream change of velocity for 50 year and 5 year floods. 
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The relation of mean velocity to discharge was plotted for 
each station as in the example of figure 1. After the discharge 
for each recurrence interval chosen had been determined from 
the flood-frequency curves, the velocity corresponding to that 
discharge was read from the velocity-discharge curves. 

The curves for recurrence intervals of 5 years and 50 years 
are reproduced in figure 4 for the stations in the upper Yellow- 
stone basin. For comparison, a few stations on the upper 
Missouri River were included in the study, and one station 
on the main stem of the Mississippi River. 

Dashed lines have been fitted by eye through the points in 
the velocity-discharge graphs of figure 4 merely to indicate 
the trend. Though there is scatter to the points, this indicates 
a remarkable conservation of the mean velocity downstream, 
and there is no definite tendency to decrease downstream in 
these rivers. A valid interpretation of these data is that in 
these streams mean velocity is essentially constant downstream. 

The graph of river slope in feet per mile plotted against 
discharge (50-year recurrence interval), which is presented for 
comparison with the velocity-discharge graphs, shows the 
marked decrease of gradient downstream from 100 feet per 
mile to about half a foot per mile at the stations studied. 
That the velocity should remain essentially constant despite 


Identification of Stream-gaging Stations Shown 
by Points in Figure 4 


. Mississippi R. at St. Louis, Mo. 
- Missouri R. at Williston, N. D. 
. Yellowstone R. nr. Sidney, Mont. 
Miles City, Mont. 
Billings, Mont. 
Corwin Springs, Mont. 
. Bighorn R. Kane, Wyo. 
Thermopolis, Wyo. 
. Greybull R. . Basin, Wyo. 
Meeteetsee, Wyo. 
. Pitchfork, Wyo. 
. Wood R. Sunshine, Wyo. 
. Clarks Fork Edgar, Mont. 
Chance, Mont. 


COWRA OK 


. Rock Cr. Joliet, Mont. 

. = nr. Red Lodge, Mont. 

. W. Fk. Rock Cr. below Basin nr. Red Lodge, Mont. 
. Red Lodge Cr. above Reservoir nr. Boyd, Mont. 

. Stillwater R. nr. Absarokee, Mont. 

. Rosebud Cr. nr. Absarokee, Mont. 


1 

12 
1 

15 
16 
17 
18 
19 
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such a decrease in slope, and under the conditions of extreme 
flood in streams of rapidly decreasing bed-material size, is of 
interest. 

Whether further study of additional data will discover 
streams in which the mean velocity rapidly decreases down- 
stream is unknown, but no such trend has been found in any of 
a large number of streams studied by the author. 


CHANGES DOWNSTREAM IN THE VELOCITY NEAR THE STREAMBED 


Though the velocity near the streambed cannot at present be 
measured directly, the theory of vertical velocity distribution 
in pipes and open channels is sufficiently well developed that 
computations of velocity near the bed are possible. Though these 
computed velocities are approximate, they allow, nevertheless, 
an indication of the downstream trends which are under dis- 
cussion here. 

The velocity distribution in the vertical tends to follow a 
logarithmic law down to a position very close to the bed itself 
where the flow is laminar, and within the laminar sublayer the 
vertical changes are linear. 

Bed velocity has been defined in a number of ways. Rubey 
(1938) considered it to be the velocity at the laminar boundary. 
Einstein (1950) discussed the movement of bed load in terms 
of the velocity in the “bed layer” which he found to be about 
2 grain diameters in thickness. 

The theory of turbulent flow derived from the study of pipe 
flow was applied to wide rivers by O’Brien (1937). Langbein 
(1940) further developed this concept and modified the O’Brien 
equations to a form which relates the velocity at a given dis- 
tance above the bed (but outside the laminar sublayer) to the 
mean velocity. As O’Brien had previously done, Langbein 
showed that actual measurements of velocity distribution gave 
results similar to those from the O’Brien equation. On the basis 
of measurement data, Langbein presents a dimensionless graph 
of the velocity distribution in the vertical which related F(z/d) 
to (z/d) where 


F(z/d) (2) 


\ 


V is the velocity at distance z from the bed, d is the depth, 
s is energy slope, V,, is the mean velocity in the vertical, and 
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F(z/d) signifies a function of (z/d). Values of the function 
F(z/d) may be read off Langbein’s graph for given values of 
z and d. For concurrent values of V,,, d, and s representing 
a given cross section, the vertical velocity distribution may be 
computed. 

Leopold and Maddock (1953) have studied the changes in 
depth, mean velocity, and water-surface slope in the downstream 
direction for numerous rivers. The downstream changes have 


been expressed as functions of discharge of a given frequency. 


From their data, which were in a form similar to those presented 
in figure 2b of the present paper, representative curves of the 
downstream changes of V,, d, and s as functions of mean 
annual discharge have been chosen and are defined by the 
mutually consistent equations below. These relations may be 
considered to describe the characteristics of a typical midwestern 
river, though no single river is exactly like the composite repre- 
sented by the equations. 


d= 0.16 Q°*° (3) 
= (4) 
(5) 


where d is depth in feet, V,, is mean velocity in feet per second, 
s is water-surface slope (approximately equal to energy slope), 
and Q is mean annual discharge in cfs. 

For a given value of Q, values of V,,, d, and s provided 
by equations 3-5 were substituted in equation 2, and using 
Langbein’s graph, the velocity V at different distances z above 
the bed was computed. The computed velocities at three dis- 
tances above the bed are shown in table 1 for several discharge 
values representing progressive downstream changes. Complete 
vertical-velocity curves for four of these assumed discharges 
are plotted in figure 5. 

The four diagrams in figure 5 represent plots of velocity 
against stream depth at four places along the length of the 
assumed composite river, The logarithmic decrease of velocity 
with depth in a given cross section is apparent from the shape 
of the curves. The Langbein empirical relation between F(z/d) 
extended only to z/d value of 0.1 and had to be extrapolated 
somewhat to compute the values of velocity near the bed for 
the deeper sections, but because his relation was a straight line 
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on logarithmic paper and parallel to the O’Brien theoretical 
relation, the extrapolation appeared to be justified. 

The measured velocities in the canal data used by Langbein 
gave the average position of the mean velocity at 0.67 the 
distance from the water surface to the streambed, whereas 
in most rivers it usually occurs at about 0.60 the depth. This 
slight difference does not materially affect the shape of the 
vertical-velocity curves. 

It can be seen in both figure 5 and table 1 that, progressing 
downstream, depth increases faster than does mean velocity. 
Moreover, the velocity at any specific distance above the bed 
also increases downstream, as can be seen in the table. In the 
case of the particular river characteristics used here as typical, 
the downstream rate of increase of velocity 0.02, 0.05 and 0.10 
foot above the bed is about the same as the downstream rate 
of increase of mean velocity. 

On figure 5 dashed lines have been drawn between the points 
representing the velocity at 0.05 feet above the bed and the 
mean velocity. The slope of the dashed line then represents 


the shear av between these two depths. The changing slope 


Ad 


| 


/ 


Mean Velocity Vm _ 


Velocity, in feet per second 


Fig. 5. Vertical velocity distribution computed for four sections repre- 
senting progressive positions downstream in a river. Note that velocity 
near the bed as well as mean velocity increases downstream. Dashed lines 
connect velocity vectors at 0.05 feet above bed and the mean velocity. 
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of these dashed lines from the upstream section to the down- 
stream section shows that this shear decreases downstream. In 
table 1, values of the shear between points 0.05 and 0.1 feet 
above the bed are listed and it can be seen that these shear 
values also decrease downstream. 

The finite values of shear indicated by the righthand column 
in table 1 are presented only to provide a better visual picture 
of the physical nature of shear stress in the vertical profile 
of velocity. They represent merely examples of shear quantities 
and it is not implied that these are particularly important 
depths through which to measure shear. 

The total shear on the bed, t, is proportional to the product 
of slope and hydraulic radius ; for wide channels where hydraulic 
radius is approximately equal to mean depth, r= ds. As indicated 
by the equations (3) and (5) for the composite river discussed 
earlier and for many of the individual rivers studied, slope 
decreases downstream slightly faster than depth increases. Thus 
the product depth x slope, to which total bed shear is propor- 
tional, decreases slightly downstream. 

It appears, then, that in numerous rivers and in a synthetic 
composite river described here, total shear on the bed and 
shear between different levels near the bed decrease in the down- 
stream direction. This decrease may obtain despite a downvalley 
increase of mean velocity and of velocity near the bed. 


RELATION TO COMPETENCE AND LOAD 


It is well known that the particle size of material composing 
the bed of a river decreases downstream owing to abrasion and 
sorting, and that there is a concomitant decrease in slope. It has 
long been supposed, therefore, that the velocity of a stream 
decreases downstream. 

The reasoning on which this supposition is based may be 
summarized as follows: The downstream decrease in caliber of 
load caused by abrasion and sorting means that a smaller veloc- 
ity is required to carry the small particles near the stream 
mouth than the large particles in headwater reaches. The smaller 
velocity required downstream is attained by progressive flatten- 
ing of the slope downstream. If, in any reach, the contributed 
load is of such large caliber that the existing velocity is in- 
sufficient to move it, deposition occurs until the slope through the 
reach of deposition is increased sufficiently to produce a velocity 
great enough to move the particles. 
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The present analysis indicates that the velocity near the bed, 
as well as the mean velocity, tends to increase downstream. It 
has been shown, furthermore, that even under conditions of 
extreme flood in steep mountain streams the mean velocity does 
not decrease downstream. It is necessary, then, to revise our 
concept of how the mutual adjustment of slope, roughness, 
velocity and caliber of load are interrelated in the downstream 
direction. 

In the present state of knowledge it is not possible to sub- 
stitute as simple and readily understandable a concept as the 
one to which we have been accustomed. A few relations which 
may be necessary and useful for further understanding can be 
stated now, but formulation of a complete picture must await 
further work. 


In streams carrying primarily coarse gravel, it is safe to 
assume that bed roughness is governed primarily by particle 
size. Downstream, where the bed is composed of fine sand, silt, 
or clay, bed roughness is materially affected by the configura- 
tion of the bed in the form of transitory riffles, waves, and dunes. 

The vertical-velocity profile is dependent in part on the bed 
roughness. Roughness affects the rate at which velocity de- 


creases as one approaches the bed—that is, the shear in the 
moving water close to the bed. This shear is intimately related 
to the production of turbulence, which in turn is responsible 
for the movement of bed particles. 

Movement of bed particles is accomplished by the energy of 
turbulent eddies immediately adjacent to the particles. A similar 
principle applies to the quarrying of large blocks of bedrock 
exposed in the stream bed and to the movement of large boulders 
or coarse gravel on the bed. On the basis of keen observation, 
Matthes (1947) described how kolk, a form of macroturbulence 
consisting of strong intermittent vortices, is a primary mecha- 
nism for plucking and for movement of large bed material. 
Vortical motion is probably, in most cases, a phenomenon aris- 
ing from shear stress. 

Therefore, a change in roughness, through its effect on shear, 
would alter the competence even without any concurrent change 
in mean velocity. Similarly, competence may be changed without 
necessarily altering slope. 

If roughness decreased downstream, shear would tend to de- 
crease downstream. This implies a downvalley decrease in 
competence. 
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The relation between shear and bed roughness and between 
these same factors and competence and load are beyond the 
scope of this paper. The importance of these relationships has 
been brought out by the many studies of bed load movement 
in which the grain size and shear factors appear prominently. 
The present paper merely calls attention to a tendency for 
certain shear functions to decrease downstream. The tendency 
for particle size to decrease downstream is well known. There 
is implied a downstream decrease of competence but any defin- 
itive discussion of this aspect requires further detailed study. 

From the foregoing analysis it is clear that we may no 
longer consider stream adjustment in terms of a downstream 
decrease in velocity. It is equally important to recognize, how- 
ever, that the introduction of the concept of decreasing shear 
downstream does not mean that we may simply substitute de- 
creasing shear or roughness for decreasing velocity in our old 
concept of stream adjustment. As can be seen in the Manning 
equation, a change of roughness, n, requires a change in the 
interrelations of velocity, depth, and slope in order that the 
equation be satisfied. Moreover, if stream adjustment is con- 
sidered in terms of the eight variables, discharge, width, depth, 
velocity, slope, roughness, load, and caliber of load, as was 
done by Leopold and Maddock (1953), then there are not one 
but eight simultaneous equations involved in the description of 
the interrelations. 
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